Distinguishing n Hamiltonians on C n by a single measurement 
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If an experimentalist wants to decide which one of n possible Hamiltonians acting on an n 
dimensional Hilbert space is present, he can conjugate the time evolution by an appropriate sequence 
of known unitary transformations in such a way that the different Hamiltonians result in mutual 
orthogonal final states. We present a general scheme providing such a sequence. 
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Controling simple quantum systems has become a large 
field of research during the last decade. Experimen- 
tal and theoretical investigations deal with the prepa- 
ration of certain quantum states, the implementation of 
unitary transformations and the design of measurement 
procedures for different quantum observables. Whereas 
the problem of optimal information gain about unknown 
quantum states is a large field of research (see e.g. jD), 
discussions about optimal estimation of unknown quan- 
tum evolutions are comparably rare. A strategy for es- 
timating an arbitrary unknown unitary transformation 
is developed in j| and in || for general quantum op- 
erations. Those approaches assume that the dynamical 
evolution to be estimated is taken from an infinite set 
of possibilities. The problem of estimating an unknown 
Hamiltonian can arise in various contextes: Assume we 
want to use a single quantum system in order to de- 
tect classical fields, e.g., a spin particle as detector for 
a magnetic field. We expose the test particle to the field 
for a certain time period and estimate the field strength 
by measuring the particle's quantum state. Assume the 
experimentalist is allowed to perform arbitrary unitary 
transformations on the test particle, expose the particle 
to the field again, repeat this several times and perform a 
single measurement at the end. What is the best proce- 
dure for estimating the field? If the set of possible values 
for the field strength is larger than n, a single measure- 
ment of the test particle can only allow estimations of 
the field. By basic quantum mechanics, it is well-known 
that a set of states is perfectly distinguishable by a single 
measurement if and only if their density matrices have 
disjoint support. Led by this simple statement concern- 
ing the distinguishability of states, we focus on the ques- 
tion of distinguishing between n possible Hamiltonians 
{H\, . . . , H n } of a quantum system on the Hilbert space 
H := C ra and show that they are always perfectly distin- 
guishable provided they do not only differ by an additive 
constant. We assume that the experimentalist is allowed 
to prepare the initial state, to perform definite unitary 
transformations interrupting the unknown natural evolu- 
tion and to perform an arbitrary measurement at the end. 
The assumption about the restricted set of possibilities is 
more natural than it might seem at first sight. Take the 
following model of a measurement interaction (compare 



||): on the joint Hilbert space of the measured system 
and the measurement apparatus we assume to have the 
Hamiltonian 



Hi 



where (Pj)j is the family of spectral projections of the 
measured observable and Hj are different self-adjoint op- 
erators moving the pointer of the measurement apparatus 
conditioned on the state of the measured system. Assume 
that we do not have any direct access to the measured 
system and that we are not able to change the interac- 
tion at all. The only way to use the interaction for a 
measurement procedure consists in initializing the mea- 
suring device, waiting (i.e. implementing e lGt ) and in- 
terrupting this evolution several times by implementing 
local unitary transformations on the measurement appa- 
ratus in order to get mutual orthogonal pointer states for 
different Hamiltonians Hj. 

Our considerations show that this is always possible (if 
Hj — tr(Hj) ^ Hi — tr(Hi) for i ^ j) and give a general 
rule for such a quantum algorithm. 

The algorithm consists of quite a large number of steps; 
since we are working in the Lie algebra instead of the 
Lie group our scheme requires arbitrarily many unitary 
transformations (close to the identity) in order to obtain 
the correct result with arbitrary reliability. We are con- 
vinced that there exist much simpler algorithms for par- 
ticular sets of n Hamiltonians. Whether or not there are 
general rules requiring only a few steps is unclear. Devel- 
oping short procedures for the general case might result 
in computationally hard word problems in the Lie group 
SU n , whereas our classical precomputation consists only 
in solving linear equations for the price of obtaining only 
approximative solutions. 

Firstly we present an example of n Hamiltonians which 
can be distinguished easily: Set Hj :— jD with D := 
diag(l, 2, . . . , n). By waiting the time t = 2tt jn we have 
implemented the unitary transformations e l ^ D / n _ Take 
the initial vector \ip) :— (1,1,..., 1) T . Then the states 
e ij2irD/n^ are or thogonal for different values of j since 
they are the discrete Fourier transforms of the canonical 
basis vectors of C™. In the rest of the paper we show that 
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the general problem can be reduced to this example. For 
doing so we start by developing some technical tools: 

By waiting the time t, we have implemented the 
transformation e lHt for the unknown Hamiltonian H G 
{Hi, . . . , H n }. We show that there is a procedure simu- 
lating e~ lHs for arbitrary s: Choose a finite subgroup Q 
of SU n acting irreducibly on Ti.. Then 

ueg 

is an operator commuting with every U S Q and is 
therefore a multiple of the identity operator by Schur's 
Lemma (this fact is used in decoupling techniques ||). 
Without loss of generality we assume every Hj to be 
traceless. Then one has YlueQ UHU' f = and hence 

J2ueg\{i} UHU] = ~ H - We obtain 

lim {Tl Ue g\sixUe iHt ' m U^) m = e~ im . 

m — >oo 

Set G := {1, Ui, . . . , Ui}. Then for large m we have ap- 
proximately an implementation of e~ lHt as follows: 

begin 

for k = 1 to to do 

for s = 1 to I do 
implement U s 
wait the time t/m 
implement C7J 

end. 

The possibility of implementing e lHt even for negative 
t is decisive for using Lie algebraic tools in the sequel: 
Let A be the Lie algebra of traceless self-adjoint opera- 
tors acting on Ti. 

By using the well-known formula 

lim ( e iH / m e iA / m e - iH / m e - iA / m ^ m2 — e i[HA] 
m— >oc 

we can design an algorithm simulating the unitary 

e -[H,A]s 

for arbitrary s£R,A€-4 with arbitrary small error. In 
the same way we conclude more generally: 

Lemma 1 Let T, Q : A — ► A be arbitrary (not necessar- 
ily linear) functions. Assume there exist for every set 
procedures for simulating the unitary transformations 

and for simulating 

e iQ{H)s 



with arbitrary small error for the unknown Hamiltonian 
H G {Hi, ... , Hn). Then there are procedures simulat- 
ing 

e i[Q {H),T{H)]s 

and 

e i[F(H),A\a 

for every A G A and every s G R with arbitrary small 
error. 

Obviously for every A € A we can find an algorithm 
performing i[H, A] =: Ad(H)(A). Hence we can find for 
every k € N an algorithm performing (Ad(H)) k (A). We 
conclude: 

Corollary Let T : A — > A be an arbitrary function. As- 
sume that for every required accuracy and every s G R 
there exists a procedure such that 

is implemented. Then Lemma 1 provides a scheme for 
implementing 

e tp(Ad(r(H))))(A) 

where p is an arbitrary real polynomial and A G A. 

Furthermore we will need the following Lie algebraic 
Lemma: 

Lemma 2 Let Hom(A, A) be the ring of R-linear maps 
on the real vector space A. Then there is no proper sub- 
ring of Hom(A, A) containing all the maps of the form 
B i— > i[A, B] with arbitrary A £ A. 

Proof Define D 1 := diag(l, -1, 0, . . . , 0), D 2 := 
(0,1,-1,0,...,0),..., D n _i := (0,..., 0,1,-1). Fur- 
thermore let Xjk for every unordered pair (J, k) with 
j < k be the matrix with 1 at the positions j, k and k, j 
and zero elsewhere. Let Yj & be defined in an analogue 
way with entries i and — i at positions j, fc and fc, j respec- 
tively. The set of these n 2 — 1 matrices forms a basis of 
the vector space A. Since all the basis vectors are unitar- 
ily equivalent (note the analogy to the Pauli-matrices), 
there always exists a map in the ring generated by maps 
of the form i[A, .] mapping one basis vector on the other. 
In order to show, that every map C G Hom(A, A) can 
be obtained by sums and concatenations of maps i[A, .] 
it is therefore sufficient to prove that a map M. can be 
generated with the following two properties: (1) The ker- 
nel of M. contains every basis vector except Di and (2) 
M{Di) is proportional to Y\$. Choose a finite subgroup 
S of SU n acting trivially on the vector space spanned by 
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the first basis vector |1) and irreducibly on its orthogonal 
complement \1} ± - The sub ring we are looking for con- 
tains clearly the map C(A) := (l/\S\)^2jUjAUj since 

UjAUj = e Ad( * B i\A) if Bj is chosen such that e iB i = Uj. 
Due to Shur's Lemma every operator in the image of C is 
a multiple of the identity on \1) ± . Since the trace on the 
subspacc ll) 1 - is invariant, the restriction of £ to If) 1 is 
given by B i — > tr(B)l. Explicitly one obtains 

C(A) = (1 - P)A(1 -P) + tr(PAP)P, 

where P is the projector onto \1) ± . Clearly C an- 
nihilates all the basis vectors except D\. Define M. 
by A4(B) := i[X\^, £{B)\ Easy calculation shows that 
A4(D\) is proportional to Y\,i. □ 

Since i[A, B] = lim s ^(,(e lAs Be~ lAs — B)/s we can con- 
clude that i[A, .] is an element of the vector space spanned 
by the maps (B \— > e tAs Be~ lAs ) se u- Hence we obtain: 

Corollary Let C £ Horn (A, A). Then there is a set 
of unitaries Ui , . . . , U rn and real numbers c\ , . . . , c m such 
ttet C{B)=Y, j < m c j U j BU). 

In order to obtain a constructive statement one can ei- 
ther take the approximative solution defined by the limit 
above or one can write Ad{A) as a finite linear combina- 
tion of maps (e Ad ( A>>s ) se m. by solving the corresponding 
equation for the eigenvalues of Ad(A). We conclude: 

Lemma 3 Let T : A — > A be arbitrary. If there is a 
scheme implementing e t:F ( H ) s for the unknown Hamilto- 
nian H £ {Hi, . . . , H n } then 

e iC(r(H))s 

for arbitrary C £ Hom(A, A) can be implemented with 
arbitrary small error by 

( Uie iFlH)c 1 ./k u 1 Uie ir(H)c a ./k u 1 u m e l ^ H ^ s ' k Ul) k 

where Uj are the unitaries and Cj are the coefficients cor- 
responding to C in the sense of the corollary to Lemma 
2 and k is large enough to keep the error small. 

Now we are able to construct our algorithm: Choose 
an operator G £ A with exactly two different eigenval- 
ues, called a and (3. Choose C £ Hom(A 1 A) in such 
a way that C{Hj) = XjG with Xj > and Aj ^ Xj. 
This is possible due to basic linear algebra. The map 
Ad(G) := i[G, .] has the eigenvalues ±i(a — (3) and 0. 
The spectrum of the map Ad(XjG) is hence given by the 
values ±Xji(a — /3),0. Choose a real polynomial p such 
that p(±Xji(a - /?)) = ±ji(a - (3) and p(0) = 0. Due 
to the functional calculus for the diagonalizable operator 
Ad(G) this implies 

p(Xj(Ad(G)))=jAd(G) 



By defining C := Ad(G)(A) for arbitrary A £ A\ {0} we 
obtain 

p(Ad(XjG))(A) = jC. (1) 

Now choose a map C £ Hom(A, A) such that 

C{C) = D2iT/n. 

We obtain C (p(Ad(C(H) ) ) (A) ) = jD2n/n. 

The classical precomputation for our algorithm can be 
sketched as follows: 

1. Choose an element G £ A with two- valued spec- 
trum and find a linear map C such that C(Hj) = 
XjG with different values Xj. 

2. Find a set of unitary transformations U\,...,Ui 
and a set of real numbers Cj such that C{B) = 
J2j c jUjBUj for every B £ A. This is possible due 
to the corollary to Lemma 2. 

3. Choose a polynomial p such that p(±Xji(a — /?)) = 
±j and p(0) = 0, if a, (i are the eigenvalues of G. 

4. Choose an arbitrary operator A £ A \ {0} and a 
map I such that C(p(Ad{G))(A)) = D2ir/n. Find 
a set of unitary operators V\ , . . . , V m and real num- 
bers dj such that C(B) = Y,j djVjBVj. 

Now we sketch the required sequence of quantum op- 
erations as follows: 

1. Prepare the initial state := (l/y/n)(l, . . . , 1) T . 

2. Call a subroutine performing the evolution e l i D27T / n 
if the Hamiltonian Hj is present. 

3. Measure in the basis defined by the discrete Fourier 
transforms of the canonical basis vectors of C". If 
the result is the j basis state then the Hamilto- 
nian Hj is present. 

The subroutine called in step (2) is recursively defined: 
The implementation of 

e ijD2w/n = e ijc(p(Ad(L(Hj)))(A)) 

is based on Lemma 2 by calling a subroutine simulating 

e ip(Ad(£(^)))(A) s 

for small s several times. The implementation of the lat- 
ter is based on the corollary to Lemma 2 by calling a 
subroutine for implementing 

e iC(H 3 )s 

several times (Lemma 3). 



3 



ACKNOWLEDGEMENTS 



Wc thank P. Wocjan for useful discussions and important corrections. This work has been supported by grants of 
the project Q-ACTA of the European Union. 



[1] 
[2] 
[3] 
[4] 
[5] 



Fuchs. Distinguishability and Accessible Information in Quantum Theory. LANL-preprint :mant-ph/9G01020 



Acin, E. Jane, and G. Vidal. Optimal estimation of quantum dynamics. LANL-preprint ^uant-ph/0012015 
D Ariano and P. Presti. Tomography of Quantum Operations. LANL-preprint qua nt-ph/001207. 
Kitaev. Quantum measurements and the abelian stabilizer problem. LANL-preprint quant-ph /95 1 1026 . 
Viola, E. Knill, and S. Lloyd. Dynamical decoupling of open quantum systems. Phys. Rev. Lett., 82:2417-2421, 1999. 



4 



